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I Abstract 

, It was shown, that the presence of the so-called X-matter with an 

• equation of state, which lies between limits of the strong and weak 

energy conditions, allows the variety of the cosmological scenarios in 
the relativistic theory of gravity. In spite of the fixed negative sign 
O , of the cosmological term in the field equations with massive graviton, 

it is possible to obtain the solutions with accelerated and complicated 
! loitering expansion of the universe. The numerical estimation of the 

universe's age agrees with the modern observational data if the upper 
limit of the graviton's mass is 10~ 71 g 



1 Introduction 



The so-called relativistic theory of gravity (RTG), which was inspired by A. A. 
Logunov with co-authors [|], ^], is the alternative to Einstein's general theory 
of relativity (GR) and differs from it in the cardinal point: the gravitation 
is considered as the tensor field of the Faraday-Maxwell's type in the flat 
Minkowski spacetime. The existence of the flat reference background metric 
in RTG allows of the unambiguity and clarity of the obtained solutions, 
which are similar to ones in GR but to have the property of the trivial 
extendibility and topological simplicity. The lasts exclude the existence of 
the singularity, black holes, wormholes, topological defects etc. The structure 
of the theory is tougher due to existence of the additional field equation on 
background metric, which is analog of the harmonic coordinate equation in 
GR and is introduced in order to exclude the spins 1 and for gravitons. 
The additional advantages of RTG are the existence of the conservation lows 
and the possibility of the canonical quantization. 
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Table 1: Cosmological parameters |^ 



cosmological parameters 


observational data 




68±6 


n K 


1.1±0.12 




0.35±0.1 




0.65±0.15 


TO 


12.7±3 


w 


< -0.6 



All experimental consequences of RTG are the same as ones in GR for 
experiments in solar system. But derived in || ||| cosmological scenario dif- 
fers essentially: RTG permits only flat global spacetime universe, critical 
density of the matter and recollapse of universe in the case of nonzero mass 
of graviton. The predictions of the universe's flatness as well as the absence 
of the cosmological singularity in the case of nonzero mass of graviton are 
the advantages of the theory. But the decelerated character of the cosmo- 
logical expansion and critical density of matter don't agree with the modern 
observational data. 

The modern estimations of the basic cosmological parameters are pre- 
sented in the Table 1 (Ho is the Hubble constant in km s _1 Mps^ 1 , &k,m,x 
are the normalized density parameters for the curvature, matter and cosmo- 
logical constant, respectively, To is the minimal estimation of the universe's 
age corresponding to age of oldest globular clusters in Gyr, w is the param- 
eter of the equation of state for the dominant energy (see below)). 

The domination of Qx provides the accelerated expansion of the universe 
because of the acceleration parameter 1 = &x — % L > (a is the scaling 
factor, H is the Hubble function, dote denotes the derivative on the time, 
the values are considered at the present moment) [||, (jj. 

In this work we demonstrate the possibility of variety of cosmological 
scenarios in RTG with the so-called "dark energy" (X-matter) term, that 
adjusts the RTG predictions with the modern astronomical data. 

2 Model 

We can not introduce the cosmological constant A with a positive sign in 
RTG aimed to describe the "repulsive" action of the vacuum without mod- 
ification of the field equations. Indeed the field equations containing A and 
terms describing the nonzero mass of graviton m have to produce a zero 
curvature of the efficient Riemannian spacetime, when the tensor of the 
energy-momentum of matter T™ = 0. That results in A = — m 2 and, as 
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consequence, "attracting" action of the vacuum (or zero cosmological con- 
stant for the massless graviton). As it was shown in J8|, there is the simple 
possibility to modification of the basic field equations due to introduction 
of the massive scalar field that results in the inflation scenario and accel- 
erated expansion of the universe. Such approach is elegant, but here we 
will consider the pure phenomenological alternative taking into considera- 
tion "dark energy" (X-matter), which does not cluster on the small scales, 
is weakly coupled to ordinary matter and has phenomenological equation of 
state px = wxpx (p is the pressure, p is the density), where wx lies between 
(dust of the "usual" matter) and -1 (true cosmological constant in GR). 
There exist the different candidates for such X-matter. The consideration of 
these models is not matter of this work (for overview see ||) and we restrict 
oneself to the phenomenological equation of state with constant negative wx- 
The Logunov's field equations are: 

G>T - ^(C + g mk ikn ~ \C9 P %k) = -8vrT n "\ (1) 



D m g mn = 0, (2) 

where G™ is the Einstein's tensor defined on the efficient Riemannian space- 
time with metrics g mn , 7 mn is the metrics of the flat background Minkowski 
spacetime, D m is the covariant derivative on the background spacetime, 
~mn _ ^jzr^gmn ^ c _ q _ ^ _ ^ rp^g choice of signs in Eq. [I] is de- 
fined by the use in the calculations of the computer algebra system Maple 6 
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Under the assumption of homogeneity and isotropy of the efficient Rie- 
mannian spacetime (the background spacetime has a Galilean metric in our 
case), its interval has the following form in the spherical coordinates [§, PL3J1 : 



ds 2 = dr 2 - a(r) 2 [dr 2 + r 2 (d6 2 + sin(fl) 2 ^ 2 )] , (3) 

where r is the proper time, a(r) is the scaling factor as function of the proper 
time. Hence from Eqs. ([| - we have the evolutional equations: 
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8 lm 2 ( 2 -l) 2 (a 2 + i) 



a I 8"< T >-i— ^ (4) 



a 4 . . . , .. 1 9 / 1 
- = --7T 3p(r) + p r " ;V 1 " - 
a 6 ova 



(5) 



Then, from the phenomenological equation of state p = wp and covariant 
conservation low V ra T^ = (here V n is the covariant derivative on Rieman- 
nian spacetime), we can obtain from Eqs. (||, ||): 



'do(t)V 
a(ty 



i - n x - n 



n 



G 



a(t) 
a(t) 2 



+ 



+ 



a(t) 1+3w x 
3 



+ 



(6) 



2a(t) 4 2a{tf 



d 2 a{t) l-n x -n R _ (l + 3wx) _ 
eft 2 2a(i) 2 2a(t) 2 + 3 <^ 1 j 

Oh „ / 1 



a(t)3 M a(t) a(t) 

where t = Hq(t—to), a(t) is the scaling factor normalized to its present value, 
&R = 1S the density parameter for the radiation (including relativistic 

neutrinos) with w R = 1/3, £lx = is the density parameter for the X- 
matter with unknown negative wx, = fjp is the density parameter for 
the gravitons. For the "usual" matter with density parameter £Im = ir£jr 
we supposed wm = 0. Also, we used the cosmic sum rule Q,m + ^r + ^x = 1 
and took into consideration that the graviton's term doesn't contribute to 
this rule. As the estimation of wx can be used the extremal values of the 
so-called weak energy condition p + p >0 and strong energy condition p + 3p 
>0 



3 Discussion 

The dynamical properties of the system (||, [?]) can be found by the search of 
the zeros of the right-hand sides of these equations. Then from the Eq. || we 
have the following polynomial equation for the minimal a m i n and maximal 
a>rnax scaling factors of the universe: 



2{n x a 3 - 3w -n G a 6 ) + 2{i-n R -n x )a 3 (8) 

+ (2VL r + 3n G )a 2 -Q G = 

The dependences of a m i n and a max on Q G are shown in Fig. 1 for fixed 
Ur « 9.1 • 1(T 5 @. As it is known [§, |, f§, 

&min 7^ for nonzero m (curve 
1). For a m i n < 10~ 4 « 1/Or (that is the approximate value of scaling factor 
at the end of the radiation epoch) the dependence of the minimal scaling 
factor on &x (and Qm) is negligible and we can write for this value: 
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Fig. 1. The logarithms of the minimal (curve 1) and maximal (curves 2-5) 
scaling factors in the dependence on the logarithm of the graviton's critical 
density. Q x = 0.65 (1-4), (5), w = -1 (2), -2/3 (3), -1/3 (4) 



But the original feature of our model is the strong dependence of a max on 
the presence and state of X-matter. If its equation of state corresponds to 
the pure A-term (w = -1), there exist the infinitely expanding solutions of 
Eqs. (g, 0) for £Ig < (curve 2, see also the asymptotically growing term 
in Eq. |8|). The transition to more "moderate" w (curves 3, 4) results in the 
recollapsing behavior of the universe with decreasing a max as result of \w\ 
decrease. The value of a max in the absence of X-matter is shown by curve 5. 
For fixed w the growth of fix from 0.5 to 0.8 (see Table 1) increases a max 
but this dependence is reduced by \w\ decrease. 

Accelerated (decelerated) character of the universe's expansion is defined 
by the sign of the polynomial, which is the modified right-hand side of Eq. 



n x a 3 - 3w (l + 3w) + 20 G a 6 + (l-tl x - M R )a 3 (10) 

+2n R a 2 - 2n G 

The roots of this polynomial give the values of the critical scaling factor a cr 
(Fig. 2) determining the transition from the accelerated (decelerated) to the 
decelerated (accelerated) expansion. 
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Fig. 2. The logarithm of the critical scaling factor in the dependence on 
the logarithm of the graviton's critical density, fix = 0.8 (1), 0.65 (2, 4, 5), 
0.5 (3), w = -1 (1 - 3), -2/3 (4), -1/3 (5). Curve 6 shows a m i n , arrows are 
directed into regions of the decelerated expansion 

It is possible to distinguish following scenarios of the universe's evolution: 

I. Deceleration and recollapse. The initially accelerated expansion changes 
into deceleration with subsequent recollapse. Such behavior takes a place for 
£la > in the case of w — -1 (right-hand side group of curves 1-3), for 
r^G lying on the right of lower curve 4 in the case of w = -2/3 and for all VLq 
in the case of w — -1/3 (curve 5). 

II. Eternal acceleration. For w = -1 there is the region of the simple 
accelerated expansion. This region increases as result of the fix increase 
(region between left-hand side and right-hand side groups of curves 1-3). 

III. Loitering expansion without recollapse. This is a complicated sce- 
nario existing for w = -1 and small £Ig (left-hand side group of curves 1-3). 
The initially accelerated expansion changes into deceleration, which changes 
into an eternal acceleration not long before the present time. The region of 
the deceleration increases as result of the £lx decrease (transition from the 
left-hand side curve 1 to 3). 

IV. Loitering expansion with recollapse. For w = -2/3 the behavior is 
the similar to III, but the long acceleration changes into deceleration and 
recollapse (curves 4). The acceleration can be very long for small Qg- 
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Fig. 3. The scenarios I - V. tt G = 0.8 (I), 0.5 (II), 10" 12 (III), 10" 3 (IV), 
0.62 (V); w = -1 (I - III, V), -2/3 (IV). fl x = 0.65. Recollapse for scenario 
IV is located outside of Fig. 



V. Marginally loitering. On the right of the region II for w — -1 there is 
the narrow region near CIq = where the accelerated expansion changes 
into decelerated one with the subsequent transition to eternal acceleration. 

The examples of the described scenarios are shown in Fig. 3. 

So, in the presence of X-matter there exist the scenarios with acceleration 
far from a m i n (II - V). As the main criterion for the selection of appropriate 
scenario we consider the demand of the existence of the radiation dominated 
epoch. The last restricts the value of CIq so as a m i n < 10~ 4 Then 
from Eq. U one can obtain Sl G < lO -11,7 that gives for the dimensional 
graviton's mass m < 10~ 71 g if Ho = 68km s -1 Alps' 1 . This estimation is 
essentially lower than in [f|]: m < 4.5 * 10 -66 g, because of the last value 
was obtained on the basis of criterion a max > 1 (i. e. the present epoch 
corresponds to expansion for VLx = 0). Obviously, that our criterion is more 
rigid. 

The obtained estimation allows to prefer the scenarios III or IV for the 
case of the accelerated expansion. The next step is the calculation of the 
universe's age, that can be made by the integration of Eq. ||: 




(11) 
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Table 2: Universe's age for scenarios III, IV 



n x 


t (w = -1), Gyr 


r (w = -2/3), Gyr 


0.5 


11.6 


11.2 


0.65 


12.9 


12.2 


0.8 


15 


13.8 



where 

_ 1 - fix - [hi _ 

a + a 1 ^™ a? 

^( a2 + i-2^)- 

It is convenient to divide the integration on the two part: radiation 
dominated epoch and matter dominated epoch. For the first interval one 
can use the approximation, which produced Eqj9|. Then the duration of 
the radiation dominated epoch is ~ 10 5 yr with slow dependence on J7g for 
dmin *C The further integration for the matter dominated epoch gives 

the results, which are presented in Table 2 (Hq = 68km s _1 Mps~ l ). The 
obtained values of the universe's age agree with the modern observational 
data (see Table 1) and don't depend on the graviton's mass. 



4 Conclusion 

The introducing of the A-term in the Lagrangian of RTG demands A = — m 2 
that contributes to the deceleration of the universe's expansion. However, 
modern observations suggest the accelerated expansion at the present epoch. 
We showed that without modification of the field equation Q the insertion 
of X-matter with the equation of the state, which is close to the extremal 
weak energy condition p + p = extends the class of the cosmological mod- 
els in framework of RTG. As result, there exist the cosmological scenarios, 
which agree with the observational data. We demonstrated, that in the case 
of m < 10~ 71 g and — 1 < w < —2/3 the minimal size of the universe, its 
age and acceleration at the present epoch don't contradict with the models 
considered in the framework of GR. But the existing data don't allow to 
choose between different scenarios because the similar present behavior re- 
sults in the essential difference in the future: eternal expansion for w = — 1 
(or m = 0), and recollapse for w > — 1. Moreover, the evolution, as rule, be- 
haves the complicated loitering character with alternation of the acceleration 
and deceleration. That demands the additional investigations in this direc- 
tion, especially, the search of RTG consequences for the universe's structure 
formation. 
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This work was carried out by the use of the computer algebra system 
Maple 6. We suppose, that the algorithmization of the RTG basics can 
stimulate the investigations in this direction. The corresponding commented 



program can be found and downloaded on 

Author is Lise Meitner fellow at the Technical University of Vienna 
(project M611). 
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